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INTRODUCTION 


For  many  non-linear  problems  it  is  generally  necessary  to  obtain 
solutions  by  either  numerical  or  approximate  methods.  In  a  recent  paper 
(1)  it  was  noted  that  typical  examples  of  these  problems  are  those  in¬ 
volving  laminar  boundary  layer  flows  and  non-linear  heat  conduction  as 
investigated  by  Biot's  variational  method.  In  that  paper  use  of  a  trans¬ 
cendental  approximation  in  laminar  boundary  layer  analysis  was  examined. 

In  the  present  paper,  use  of  the  same  transcendental  approximation  in 
non-linear  heat  conduction  is  discussed. 

In  previous  studies  of  typical  one-dimensional,  transient  heat 
conduction  problems  with  non-linear  boundary  conditions  (2,  3)  the 
temperature  distribution  has  been  represented  subject  to  a  two-fold 
approximation.  First,  the  boundary  condition  infinitely  distant  from  the 
surface  is  brought  to  a  finite  distance  from  the  surface,  this  distance 
being  called  the  penetration  depth;  then  the  temperature  is  approximated 
by  a  polynomial  in  the  region  from  the  surface  to  the  penetration  depth, 
being  regarded  as  constant  at  greater  depths.  For  the  cases  considered 
previously ,  it  is  possible  to  obtain  quickly  at  least  the  asymptotic 
solutions  for  short  and  for  long  times.  However,  it  is  not  clear  that 
the  solutions  will  be  roughly  of  equal  accuracy,  since  the  polynomial 
profile  used  may  be  less  suited  to  represent  the  true  profile  at  one  of 
the  (asymptotic)  times  than  at  the  other.  The  profile  function  used 
here  is  adaptable  to  the  extent  that  the  profile  shapes  obtained  for  each 
asymptotic  case  need  not  be  the  same.  Further,  it  is  known  that  the  actual 
profile  must  be  some  form  of  transcendental  function,  and  a  transcendental 
approximation  can  have  a  closeness  of  indefinitely  high  order  to  the 
exact  solution. 

In  this  paper  a  particular  example  of  a  transient  conduction  problem 
with  a  non-linear  boundary  condition  is  considered,  for  which  other  solu¬ 
tions  have  been  obtained.  This  example  consists  of  unsteady  one-dimensional 
conduction  in  a  semi-infinite  slab  with  the  boundary  condition  that  the  heat 
flux  at  the  surface  is  proportional  to  the  n-th  power  of  the  surface  temper¬ 
ature  and  with  the  initial  condition  that  the  slab  temperature  is  uniform, 
forming  the  reference  zero,  i.e.  the  zero  of  an  appropriate  empirical 
temperature  scale. 


Solutions  for  this  form  of  boundary  condition  are  useful  for  a 
variety  of  problems.  One  practical  problem  for  which  they  have  been 
used  is  transient  conduction  due  to  unsteady  radiation  in  an  enclosure. 
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NOTATION 


A 

a 

B 

b 

C 

c 

D 

-Ei(-u) 

Ep(u) 

F 

f 

gk 

H 

I 

.m  - 
1  erfc  n 

J 

K 

k 

M 

M 

m 

m 

n 

n 

P 


S 

T 

t 

u 

u 

o 

V 

V 


coefficient 

profile  parameter,  function  of  n 
coefficient 

profile  parameter,  function  of  n 

coefficient 

specific  heat 

dissipation  function;  also,  a  coefficient 

exponential  integral 

general  exponential  integral 

heat  flux  at  x  =  0 

coefficient 

constants 

heat  flow  field;  H,  heat  flux  field 
a  specific  integral 

m-th  repeated  integral  of  the  error  function 

coefficient 

coefficient 

thermal  conductivity;  also,  running  index 
function  of  u 

o 

coefficients 
running  index 
exponent 

unit  normal  vector 
running  index 
thermal  force 

i-th  generalized  coordinate;  q^,  surface  temperature; 

q^,  penetration  depth 
surface 

arbitrary  constant  temperature 
time 

exp  (a  +  bn) 
exp  a 

a  thermodynamic  potential 
volume 
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x»  y>  z 

pro  n 


“o’  eo’  Yo 

r(l+n) 


Y 

e 

C 

n 

0 


p 

o 


T 


X 

♦ 


space  coordinates 

dimensionless  temperature  profile  approximation 

coefficients 

gamma  function 

Euler* s  constant 

an  arbitrary  small  number 

dimensionless  space  variable 

general  variable 

temperature 

density 

scaling  factor 
dimensionless  time 
dimensionless  penetration  depth 
dimensionless  surface  temperature,  x  =  0 
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Biotas  Variational  Method  for  Transient  Conduction 

A  variational  principle  for  analysis  of  heat  conduction  has  been 
described  extensively  by  Biot  (t,  5),  who  has  also  given  some  examples  of 
its  use.  An  alternative,  complementary  method  for  deriving  the  principle 
was  described  by  Boley  S  Weiner  (6).  In  essence,  it  is  supposed  that  one 
can  consider  a  heat  flow  vector  field  H  to  exist  within  a  body,  such  that 
the  time-rate  of  change  H  is  the  heat  flux  across  an  area  normal  to  H  . 
From  the  First  Law,  in  the  absence  of  work,  one  has  that 


cp0  =  -div  H 

[1] 

In  particular,  if  the  heat  flow  field  can  be  expressed  as  a 
n  "generalized  coordinates"  q^  (t)  ,  so  that 

function  of 

H  =  H  (q.;  x,  y,  z,  t) 

C2] 

then  it  is  possible  to  write  a  variational  principle  in  the 

form 

9qi  H,  1 

[3] 

where 

Qi  =  0  n*(3H/Bqi)  dS 

[4] 

and,  common  to  the  n  equations  [3], 

V  =  J  (l/2)cp  e2  dv 

C5] 

”  ■  j  k  "2  dv  • 

[6] 

This  variational  principle  is  equivalent  to  the  heat  conduction  equation 
in  an  isotropic  medium.  When  the  temperature  field  is  one -dimensional 
in  space,  it  and  the  heat  flow  field  are  related  simply.  In  order  to 
utilize  the  method,  some  particular  form  of  temperature  profile  is  assumed, 
and  a  sufficient  number  of  generalized  coordinates  assigned  to  describe 
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it.  From  eqn.  [1]  the  heat  flow  field  corresponding  to  the  assumed 
profile  can  be  determined,  and  subsequently  the  dissipation  function  D 
and  the  potential  V  . 

As  an  alternative  to  determining  the  n  generalized  coordinates 
from  eqn.  [3],  it  is  possible  to  determine  one  through  use  of  a  compatibi¬ 
lity  condition  on  the  surface  heat  flux  (e.g.,  by  forcing  Hq  for  a  semi¬ 
infinite  body  to  satisfy  exactly  the  heat  flux  boundary  condition), 
leaving  (n  -  1)  coordinates  to  be  determined  from  the  variational 
principle. 


The  Profile  Function 


The  function  for  profile  approximation  introduced  by  Richardson 
(7)  is 


pro  x]  =  exp  [-  exp  (a  +  bn)]  /  exp  [-  exp  a]  [7] 

and  this  function  satisfies  the  conditions 

pro  n  =  l  at  n  =  0 

pro  n  **  0  as  n  **  00  . 

This  profile  function  is  a  two-parameter  profile,  and  in  general  two  in¬ 
dependent  equations  are  required  in  order  to  determine  the  relevant 
parameters. 


The  Differential  Equations  for  the  Profile  Function 

The  equations  found  by  Biot’s  method  depend  upon  the  specific 
temperature  profile  which  has  been  assumed.  It  is  convenient  to  write 

n  =  x/q2  , 
u  =  exp  (a  +  bn)  , 
and  uq  =  exp  a  , 
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so  that  the  profile  assumed  is 


8  =  qx  exp  (~u)  /  exp  (~uq)  .  [8] 

It  is  also  assumed  that  the  heat  flux  boundary  condition  F  belongs  to 
the  class  of  functions 


F  =  f0on  =  f  qxn  •  [9] 

It  is  convenient  to  form  dimensionless  groups  for  q^  and  t  ,  but  it 

is  found  that  the  parameters  of  the  problem  do  not  provide  enough  quanti¬ 
ties  and  it  is  necessary  to  introduce  an  arbitrary  temperature  or  length 
(which  subsequently  cancels  out  in  the  solutions).  If  the  former  is  chosen 
and  denoted  by  ,T  ,  the  following  dimensionless  groups  can  be  established: 


*  a  qx  /  t 

x  =  q2f  Tn_1  /  k 
T  =  t  f2  T2(n_1)  /  kcp 


[10] 


The  heat  flow  field  is  given  by 


H  = 


cpqnq9 

r - M — r-  Ei  (-u)  , 

b  exp  C-u  )  ’ 


[11] 


so  that  by  determining  eqn.  [3]  for  qx  it  can  be  shown  that 


J  tp  +  ( 3/2  )Ai^xX  +  Ai|>x  =  0 


[12] 


while  the  surface  flux  compatibility  condition  gives  that 


'f'X  +  • 


[13] 


These  equations  have  the  same  form  as  those  obtained  by  Richardson  (3), 
but  the  coefficients  differ  because  of  the  different  assumed  profile. 
The  coefficients  which  arise  here  are 
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[14] 

[15] 

[16] 

[17] 


In  deriving  eqn.  [12],  advantage  is  taken  of  the  fact  that  B  =  A/2  , 
which  can  be  shown  by  noting  that  n  =  (In  u  -  a)/b  ,  transforming  to 
u  as  the  variable  of  integration,  and  integrating  by  parts. 

Equations  [12]  and  [13]  can  be  solved  easily  for  the  cases  n  =  0 
(for  which  there  is  also  an  exact  solution)  and  0  <  n  <  1  .  The  solution 
for  n  =  0  (i.e.  constant  heat  flux)  is 


X  =  Yct 


1/2  _ 


Yo  (2  -  n)  1/(n  1)  1/2(1  -  n) 

*  =  IT  2(1  “nT  T 


It  is  to  be  noted  that  this  solution  cannot  be  used  near  n  =  1  .  (In 
(3)  these  sets  of  solutions  were  incorrectly  indicated  to  be  asymptotic 
solutions  instead  of  complete  solutions,) 


8 


Evaluation  of  Coefficients 


The  next  step  in  completing  the  solutions  is  to  evaluate  the  co¬ 
efficients  A,  J  and  K  •  These  coefficients  can  be  obtained  as  functions 
of  a  and  b  .  Finally  a  criterion  can  be  established  for  obtaining  a 

and  b  ,  thereby  completing  the  solutions. 

3 

It  may  be  noted  that  Ab  ,  Jb  and  K/b  are  each  functions  of  uq 

only 9  and  therefore  of  a  only.  It  can  be  seen  from  their  definitions 

that  Jb  and  K/b  can  be  obtained  from  tables  of  exponential  integrals 

3 

using  straightforward  manipulations.  The  integral  for  Ab  cannot  be 
evaluated  directly,  but  the  integral  can  be  expanded  in  an  infinite 
series  which  can  be  integrated  term  by  term.  It  was  found  that  the 
solutions  required  use  of  values  of  uq  roughly  in  the  range  of  3.0  to 
8,0  .  For  small  values  of  the  argument  the  exponential  integral  can  be 
written 


Ei(-u)  =  y  +  ln(u)  -  u  +  (l/2!2)u2  -  (l/3!3)u3  +  ...  [20] 

This  series  converges,  but  for  values  of  u  >  1  the  terms  initially 
increase  in  magnitude,  making  it  necessary  to  take  many  terms  to  obtain 
accurate  results.  For  large  values  of  u 

Ei(-u)  =  e  U  {u~X  -  u  2  +  2u  ^  -  (3!)u  4  [21] 


This  second  series  is  asymptotic,  the  terms  initially  converging  and  then 
diverging. 

It  can  be  seen  that 


A  = 


exp  (2u  )  r«> 
o 


u 


(Ei(-u) } 


2  du 
u 


exp  (2u  ) 
o 


I  ,  say. 


[22] 


For  values  of  uq  4  the  asymptotic  expansion  [21]  was  substituted, 
giving 
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■I 
-I 


-2u  f  y  /  _  v iri'f'  1  ,  ,vt  -m>2  du 
e  {  l  (-1)  (m  -  l)!u  }  — 

u  m=l 

o 


OO  00 

-2u  v  t  i  -k  , 

e  l  (-D  gw  u  du 

u  k=3 

o 

00  / 

l  (-l)k  g.  U_k  exp  (-2u)  du  , 

=■*  K  J„ 


[23] 


k=3 


where  values  of  are  g3  =  1,  g4  =  -2,  gg  =  5,  gg  =  -16,  g 

Now,  the  general  exponential  integral  is  defined  (8)  as 


=  64, 


g8  =  -312,  etc 


E  (u)  = 
P 


x  ^  exp  (-xu)dx  , 


[24] 


so  that  the  typical  integral  term  of  I  is 
|  u  k  exp  (-2u)du  =  u^  k  E^(2uq) 

J  11 


whence 


I 


-u 

o 


w  k 

l  (-1)  g 

k=3 


k 


(2u  ) 
o 


so  that  the  coefficient  A  could  be  evaluated  as  a  series  of  general 
exponential  integrals 


u  exp  (2u  )  “  k  _k 

A  =  -  3  l  (-1)"  gk  \  Ek  <2uo)  . 

jd  K**  o 


[25] 


In  computation  the  series  was  summed  to  the  smallest  term;  if  the  consequent 
maximum  error  was  considered  too  large,  Euler's  transformation  was  applied, 
and  further  terms  taken  until  the  maximum  error  was  reduced  to  a  satisfac¬ 
tory  magnitude. 


Determination  of  Constants  for  the  Profile  Function 

The  coefficients  A,  J  and  K  of  the  differential  equations  were 

determined  as  functions  of  u  and  b  in  the  last  section.  The  values 

o 
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of  uq  can  be  found  by  application  of  the  heat  flux  boundary  condition, 
eqn.  [9]: 


F  =  f  6n  =  f  q”  =  -k(30/3x)  . 

0.1  o 

Differentiation  of  eqn.  [8]  gives 
F  =  k  qx  b  Uq  /  q2  , 
which  leads  to 

^n’1  X  /  b  uq  =  1  . 

Substitution  of  the  solutions  [19]  gives 

bu  A  K  /  J  =M=  2(2-  n)/(3n  -  5)  .  [26] 

o 

The  left  hand  side  of  eqn.  [26]  is  a  function  of  uq  only,  but  the  value 
of  uq  corresponding  to  a  specific  value  of  n  cannot  be  extracted 
directly.  Values  of  M  were  computed  for  values  of  u=  3.0  (1.0)  8.0 
with  use  of  tables  and  the  asymptotic  expansion,  eqn.  [21],  and  for  values 
of  uq  £  4.0  by  integration  using  eqn.  [20].  In  calculation  at  uq  =  4  , 
one  of  the  series  used  in  determining  A  had  to  be  taken  to  twenty-five 
terms  to  evaluate  the  integral  to  eight  significant  digits.  Checks  were 
made  by  overlapping  and  use  of  alternative  series.  In  order  to  determine 

values  of  u  for  specified  n  ,  it  was  found  convenient  and  adequate  to 

.°  4 

use  linear  interpolation  between  values  of  M  vs.  log  uq  .  It  was 

estimated  that  interpolation  errors  would  correspond  to  less  than  +_  0.001  M  , 

Once  values  of  u  have  been  found,  values  of  a  can  be  determined  immedi- 

ately.  The  values  of  b  remain  to  be  found. 

For  the  polynomial  approximation  of  the  temperature  distribution  the 

penetration  depth  is  defined  as  the  point  where  0=0, i.e.  x  =  q2  . 

With  the  profile  function  0  0  as  x  ■+  *>  and  hence  q2  must  be  defined 

differently.  It  is  convenient  to  put  x  =  where  0/q1  =  e  ,  where  e 

is  an  arbitrary,  very  small  number.  From  the  profile  function  it  follows 

that 
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b  =  ln(l  -  e“a  In  e)  . 


[27] 


It  was  decided  to  use  z  =  0,01  here.  It  must  be  emphasized  that  the  value 
of  b  depends  upon  the  value  of  z  .  With  this  convention,  it  is  not 
significant  to  compare  directly  the  coefficients  of  the  penetration  depth 
solutions  from  use  of  polynomial  and  transcendental  profiles. 


Results  and  Accuracy  of  Computations 


Constant -flux  solution  The  values  found  for  the  constant-flux 

solution  were  u  =  2.95  , 
o 


=  1.130  x 


1/2 


[28] 


X  =  3.18  t 


1/2 


which  can  be  compared  with  Richardson’s  (3)  parabolic  approximation 


4  =  1.157  t1/2  ;  X  =  2.59  i1/2 


and  with  the  exact  solution 

1/2 

=  1.1284  x  • 

1  /O 

It  can  be  seen  that  the  coefficient  for  t|//T  using  the  profile  function 
differs  from  the  exact  solution  by  about  0.2  per  cent  and  from  the  para- 
bolic  approximation  by  about  2.4  per  cent. 

Solutions  for  n  >  0  A  set  of  solutions  was  computed  for  a  range  of 
n  where  the  computations  did  not  become  inordinately  long.  In  Table  I 
values  of  a  and  b  are  provided,  and  in  Table  II  values  of  the  surface 
temperature  coefficients  (parabolic  profile),  ty/x 

(transcendental  profile)  and  penetration  depth  coefficient  (transcendental 
profile)  are  listed,  together  with  the  percentage  differences  of  the  sur¬ 
face  temperature  coefficients. 

For  the  sake  of  comparison,  two  temperature  profiles  have  been  plotted 
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in  Fig.  1.  The  scale  of  slab  depth  n  has  been  adjusted  by  using  scaling 
factors  a  so  that  both  profiles  have  the  same  slope  at  n  “  0  •  It  can 
be  seen  that  the  profile  for  the  constant-flux  solution  has  a  somewhat 
different  shape  from  that  for  the  solution  with  n  =  0.7  .  Profiles  for 
intermediate  values  of  n  fall  between  those  shown. 

Another  comparison  has  been  made  in  Fig.  2  for  the  constant-flux 
solution.  The  ordinate  is  of  temperature,  normalized  such  that  the 
temperature  on  the  axis  (which  corresponds  to  the  body  surface)  is  always 
unity;  the  abscissa,  the  slab  depth,  has  been  adjusted  such  that  all 
temperature  profiles  shown  have  a  slope  of  unity  at  the  body  surface.  The 
middle  solid  line  represents  the  profile  of  the  analytic  solution  of  the 
constant  flux  case;  two  other  solid  lines  represent  the  parabolic  and  the 
simple  exponential  approximations  respectively.  The  dashed  line  represents 
the  profile  approximation  used  here.  Since  this  lies  so  close  to  (and 
crosses)  the  analytic  solution,  only  parts  can  be  shown  without  obscuring 
the  analytic  solution.  This  figure  demonstrates  well  that  the  profile 
function  used  here  is  a  closer  approximation  than  are  the  others  cited. 

The  computations  using  the  exponential  integrals  were  based  upon  the 
tables  of  Pagurova  (8),  with  values  for  integrals  beyond  the  range  of  the 
tables  generated  using  relations  given  in  the  same  reference.  Computations 
for  other  quantities  described  here  were  made  using  the  handbook  of 
Abramowich  and  Stegun  (9),  with  checks  being  made  in  relevant  references 
given  therein.  The  accuracy  of  the  computations  varies  with  n  in  the 
solutions,  but  it  is  believed  to  be  not  worse  than  0.1  per  cent  in  any 
case. 


Analytic  Solutions  for  n  =  1.0 

The  analytic  solution  for  n  =  0  ,  i.e.  constant  heat  flux,  is  well 
known.  The  problem  can  also  have  an  analytic  solution  at  n  =  1.0  if 
the  initial  condition  described  in  the  Introduction  is  slightly  altered 
so  that  the  slab  temperature  is  uniform  but  above  the  reference  zero  by  an 
arbitrary  amount.  Without  this  alteration  the  solution  recovered  is  the 
trivial  case  0(x,  t)  =  0  .  The  arbitrary  uniform  temperature  can  be 
chosen  as  the  unit  temperature  of  the  0  temperature  scale;  the  solution 
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is  proportional  to  it.  The  solution  can  be  written 


e(x,  t)  - 


i  = 


00 


l 


m+1 


^(m+l)/2  i^^erfcCc/i1^2) 
im+^erfc(0) 


[29] 


with 


J2_ 

Jn 


r(J+  l) 


m+l 


r<2ii  t  l) 


T  = 


tf 

kcp  * 


C  = 


xf 

2k 


(where  f  has  been  normalized  with  respect  to  the  initial  temperature). 

The  function  im  erfcn  is  the  m-th  repeated  integral  of  the  error  function, 
for  which  useful  tables,  recurrence  relations,  expansions  and  so  forth 
exist.  At  small  times,  the  term  for  m  =  0  is  dominant  and  this  coincides 
with  the  constant  flux  solution  discussed  previously ,  For  longer  times,  the 
series  provides  a  shifting  average  of  the  repeated  integrals  of  the  error 
function.  Successive  integrals  have  profile  shapes  which  move  from  the 
central  profile  of  Fig.  2  towards  the  simple  exponential.  The  approximate 
profiles  found  here  for  0  <  n  <  1.0  have  a  trend  in  the  same  direction. 

For  very  long  times,  the  asymptotic  behavior  of  im  erfcn  as  m  00  is 
important.  This  can  be  determined  from  the  series  expansion 


im  erfcn  =  £ 

p=0 


(-l)p  np 

2m'p  P!  r(i  +  2!=E) 


in  which  the  terms  corresponding  to  p  =  m+2,  m+4,  m+6,  ...  are  under¬ 
stood  to  be  zero.  As  m  tends  to  infinity,  it  can  be  shown  that  the 
normalized  repeated  integral  of  the  error  function  (im  erfcn/im  erfcO) 
is  given  asymptotically  by  a  simple  negative  exponential  function  of  n  • 


Discussion 


It  can  be  seen  from  Fig.  1  that  the  solution  of  the  problem  does 
take  advantage,  so  to  speak,  of  the  flexibility  of  the  transcendental 
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profile  by  finding  that  different  shapes  are  appropriate  to  different 
conditions.  This  adaptive  feature  of  the  transcendental  profile  should 
assist  in  providing  results  which  are  more  accurate  than  with  fixed 
simple  profiles.  This  hope  is  borne  out  clearly  in  the  constant-flux 
solution,  where  the  accuracy  of  the  surface  temperature  coefficient  is 
improved  by  a  tenfold  order  of  magnitude  compared  with  the  simple  poly¬ 
nomial  and  exponential  profiles  listed  in  Table  I  of  LardnerTs  paper  (2). 

It  is  also  demonstrated  by  Fig.  2. 

For  0  <  n  <  1.0  no  exact  solutions  exist  with  which  the  results 
can  be  compared.  The  calculations  become  increasingly  difficult  as  n 
approaches  unity.  This  is  due  partly  to  the  need  to  generate  exponential 
integrals  beyond  the  range  of  available  tables,  and  partly  to  the  behavior 
of  the  profile  function.  In  the  profile  function  the  parameter  a  essen¬ 
tially  specifies  the  shape  of  the  profile,  while  the  parameter  b 
specifies  the  extent  of  the  profile.  Thus,  two  profiles  drawn  with  the 
same  values  of  a  and  different  values  of  b  have  the  same  shape,  but 
not  vice  versa.  As  a  -*  *>  ,  prori  exp  [-(b  exp  a)n  ].  Over  the  major 
range  of  pron  ,  this  limit  is  approached  rapidly.  Even  with  a  =  5.0  , 
pron  is  close  to  its  limit.  This  means  that  if  an  attempt  is  made  to  fit 
the  profile  function  to  a  function  which  is  close  to  the  simple  exponential, 
the  determination  of  a  becomes  extremely  insensitive;  large  changes  in  a 
produce  small  changes  in  the  profile  shape.  It  is  increasingly  difficult 
to  obtain  values  of  a  to  a  specified  number  of  significant  digits. 

It  is  noteworthy  that  the  difference  between  the  surface  temperature 
coefficients  listed  in  Table  II  increases  smoothly  from  the  value  of  2.4% 
at  n  =  0  to  12%  at  n  =  0,73  ,  the  difference  increasing  roughly  expo¬ 
nentially  with  n  .  The  difference  between  surface  temperature  coefficients 
is  sufficiently  large  to  be  significant  in  application. 

These  differences  occur  because  the  transcendental  profile  changes 
its  shape  with  n  ,  which  the  parabolic  profile  cannot  do.  It  is  very 
probable  that  the  solutions  with  the  transcendental  profile  are  more 
accurate,  but  at  present  this  cannot  be  demonstrated  directly.  However, 
it  is  possible  to  make  comparisons  with  the  analytical  solutions  for  n  =  0 
and  n  =  1.0  and  show  that  the  temperature  profile  shape  generated  here 
for  0  <  n  <  1.0  varies  smoothly  between  the  analytic  limits.  The  signi¬ 
ficant  feature  of  the  results  upon  which  a  comparison  can  be  based  is  the 
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parameter  a  ,  which  determines  the  shape  of  the  profile.  It  has  been 
noted  above  that  the  asymptotic  behavior  of  the  analytic  solution  for 
n  =  1.0  is  that  the  profile  becomes  a  simple  exponential;  in  approximat¬ 
ing  this  shape  a  becomes  large.  For  n  =  0  ,  the  short-time  solution 
is  valid  at  all  times,  so  that  in  this  instance  it  is  also  the  long-time 
solution.  The  values  of  a  for  solutions  in  the  range  0  <  n  <  1.0  are 
shown  in  Fig.  3.  In  this  it  can  be  seen  that  the  values  found  here  fit 
well  between  the  analytic  limits.  This  feature  is  absent  in  all  fixed 
profiles,  such  as  polynomials.  The  comparison  provides  conducive  evidence 
that  the  solutions  presented  here  are  considerably  more  accurate  than  the 
corresponding .solutions  available  previously. 


Summary 


(1)  Biot!s  variational  method  is  applied  to  a  problem  of  transient 
heat  conduction  in  a  semi-infinite  slab  subject  to  a  non-linear 
boundary  condition. 

(2)  A  two-parameter  transcendental  approximation  is  used  for  the 
temperature  profile.  This  approximation  has  the  advantage  that  its 
shape  is  not  fixed,  so  that  the  profile  determined  for  each  case 
considered  can  have  the  shape  most  appropriate  to  it. 

(3)  Computations  for  solutions  utilizing  the  general  exponential 
integral  are  described,  and  two  exact  solutions  are  mentioned  for 
cases  which  bound  the  examples  computed. 

(4)  Comparison  of  the  variational  solution  using  the  transcendental 
approximation  with  the  exact  solution  for  the  limiting  case  of 
constant  heat  flux  demonstrates  an  error  of  less  than  0.2  per  cent 
in  the  surface  temperature  coefficient  and  very  close  representation 
of  the  true  temperature  profile .  This  corresponds  to  a  very  consider¬ 
able  improvement  over  solutions  obtained  previously  with  other 
profile  approximations. 

(5)  Comparison  of  the  variational  solution  for  other  examples  shows 
that  the  computed  profile  shapes  have  a  uniform  variation  which  is 
in  the  correct  direction  to  merge  with  the  limiting  case  of  surface 
heat  flux  directly  proportional  to  surface  temperature. 
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(6) 


It  is  concluded  that  the  solutions  obtained  here  using  the 
transcendental  approximation  are  considerably  more  accurate  than  those 
previously  available. 
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Table  I 


Solutions  for  0  <  n  <  1 


n 

a 

b 

0.1 

1.13287 

0.90959 

0.2 

1.21696 

0.86023 

0.3 

1.31709 

0.80370 

0.4 

1.43553 

0.73998 

o 

• 

tn 

1.57704 

0.66855 

0.6 

1.75592 

0.58528 

0.7 

1.98445 

0.49042 

0.73 

2.06623 

0.45951 
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Table  II 


n 

Surface  temperature  coefficients 
#/Tl/2(l-n) 

Penetration  depth 

/  i/2 

X/T 

Difference  in 

surface  tempera- 

parabolic 

transcendental 

transcendental 

ture  coefficients 

profile 

profile 

profile 

(per  cent) 

0.1 

1.1347 

1.1053 

3.09136 

2.7 

0.2 

1.0996 

1.0658 

3.05915 

3.2 

0.3 

1.0453 

1.0064 

3.01486 

3.9 

0.4 

0.96170 

0.91769 

2.96901 

4.8 

0.5 

0.83340 

0.78593 

2.92637 

6.0 

0.6 

0.64006 

0.59309 

2.76203 

7.9 

0.7 

0.37095 

0.33493 

2.57070 

10.7 

0.73 

0.28085 

0.25068 

2.49858 

12.0 
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Figure  1.-  Comparison  of  the  shapes  of  temperature  profiles.  The  profiles  have 
been  drawn  with  different  scale  factors  to  make  the  initial  slopes  identical. 
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Figure  2.-  Comparison  of  the  exact  temperature  profile  for  the  constant  flux 
solutions  with  various  approximations.  The  approximation  function  used 
here  (dashed  curve)  gives  a  very  close  approximation  to  the  exact  profile 


Figure  3.-  Comparison  of  the  profile  shape  parameter  for  solutions  as  a  func¬ 
tion  of  n.  At  n  =  1  the  exact  solution  corresponds  to  a  -» oo^  but  even 
log  a  =  0.70  gives  a  very  close  approximation. 
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